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Abstract 

We consider the perturbation of parabolic operators of the form dt + P{x, D) 
by large-amphtude highly oscillatory spatially dependent potentials modeled as 
Gaussian random fields. The amplitude of the potential is chosen so that the 
solution to the random equation is affected by the randomness at the leading 
order. We show that, when the dimension is smaller than the order of the elliptic 
pseudo-differential operator P{x,D), the perturbed parabolic equation admits a 
solution given by a Duhamel expansion. Moreover, as the correlation length of 
the potential vanishes, we show that the latter solution converges in distribution 
to the solution of a stochastic parabolic equation with a multiplicative term that 
should be interpreted in the Stratonovich sense. The theory of mild solutions for 
such stochastic partial differential equations is developed. 

The behavior described above should be contrasted to the case of dimensions 
that are larger than or equal to the order of the elliptic pseudo-differential operator 
P{x, D). In the latter case, the solution to the random equation converges strongly 
to the solution of a homogenized (deterministic) parabolic equation as is shown 
in the companion paper [2]. The stochastic model is therefore valid only for 
sufficiently small space dimensions in this class of parabolic problems. 

keywords: Partial differential equations with random coefficients. Stochastic partial 
differential equations, Gaussian potential, iterated Stratonovich integral, Wiener-Ito 
chaos expansion 

AMS: 35R60, 60H15, 35K15. 



1 Introduction 

We consider the parabolic equation 



u^iO,x) = Uo{x), 



*Dcpartment of Applied Physics and Applied Mathematics, Columbia University, New York NY, 
10027; gb2030@columbia.edu 

1 



where P{x, D) is an elliptic pseudo-differential operator with principal symbol of order 
m > d and x E W^. The initial condition Uo{x) is assumed to belong to L^(M'^) fl L^(M"'). 
We assume that q{x) is a mean zero, Gaussian, stationary field defined on a probability 
space {fl,J^,F) with integrable correlation function R{x) = K{q{0)q{x)}. 

The main objective of this paper is to construct a solution to the above equation in 
L^(r2 X M*^) uniformly in time on bounded intervals (see Theorem [3] below) and to show 
that the solution converges in distribution as e — > to the unique mild solution of the 
following stochastic partial differential equation (SPDE) 

— + P{x,D)u-auoW = ^2) 

u(0, x) = Uq{x), 

where W denotes spatial white noise, o denotes the Stratonovich product and sigma is 
defined as 

:= (27r)'^^(0) = [ E{q{Q)q{x)}dx. (3) 

We denote by G(t, x; y) the Green's function associated to the above unperturbed 
operator. In other words, G{t,x;y) is the distribution kernel of the operator e~^^^^'^\ 
Our main assumptions on the unperturbed problem are that G{t,x]y) = G{t,y]x) is 
continuous and satisfies the following regularity conditions: 

sup / \G{t,x;y)\dx + sup t"^ / \G{t,x]y)\'^dx + sup t'^\G{t,x;y)\ < oo. (4) 

t,y Jr^ t,y J^d t,x,y 

Note that the bound is a consequence of the and L°° bounds. Such regularity 
assumptions may be verified e.g. for parabolic equations with m = 2 and d = 1 or more 
generally for equations with m = 2n an even number and rf < m. The convergence 
of the random solution to the solution of the SPDE is obtained under the additional 
continuity constraint 

sup s"' I \G{s,x,C,) — G{s,x + y^Q\dx—^Q as y — for 7 = 2(1 j. (5) 

Such a constraint may also be verified for Green's functions of parabolic equations with 
m = 2n and < m; see lemma 14.11 below. 

We look for mild solutions of ([2]), which we recast as 

u{t,x) = e-*^^^\o{x)+ [ [ G{t- s,x]y)u{s,y)o(7dW{y)ds. (6) 

Jo Jr'' 

Here, dW is the standard Wiener measure on M'^ and o means that the integral is 
defined as a (anticipative) Stratonovich integral. In section [2l we define the Stratonovich 
integral for an appropriate class of random variables and construct a solution to the 
above equation in L^(fi x M.'^) uniformly in time on bounded intervals by the method of 
Duhamel expansion; see Theorem [1] below. In section [3l we show that the solution to 
the above equation is unique in an adapted functional setting. The convergence of the 
solution u^{t) to its limit u{t) is addressed in section HI see Theorem H] below. 
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The analysis of stochastic partial differential equations of the form ([2]) with m = 2 
and with the Stratonovich product replaced by an Ito (Skorohod) product or a Wick 
product and the white noise in space replaced by a white noise in space time is well 
developed; we refer the reader to e.g. [5], El [HI [HI [13 [21] • The case of space white 
noise with Ito product is analyzed in e.g. [1]. One of the salient features obtained in 
these references is that solutions to stochastic equations of the form ([2]) are found to 
be square-integrable for sufficiently small spatial dimensions d and to be elements in 
larger distributional spaces for larger spatial dimensions; see in particular [5] for sharp 
criteria on the existence of locally mean square random processes solution to stochastic 
equations. This brings into question the justification of stochastic models of the form 

The theory presented in this paper shows that the solution to ([2]) may indeed be 
seen as the e — limit of solutions to a parabolic equation ([1]) with highly oscillatory 
coefficient when the spatial dimension is sufficiently small. In larger spatial dimensions, 
the behavior observed in [2] is different. The solution to ([T]) with a properly scaled 
potential (of amplitude proportional to for m < d) converges to the deterministic 
solution of a homogenized equation, at least for sufficiently small times. The solution 
to a stochastic model no longer represents the asymptotic behavior of the solution to 
an equation of the form ([T]) with highly oscillatory random coefficients. 

The analysis of equations with highly oscillatory random coefficients of the form ([T]) 
has also been performed in other similar contexts. We refer the reader to [T7] for a 
recent analysis of the case m = 2 and d = 1 with much more general potentials than the 
Gaussian potentials considered in this paper. When the potential has smaller amplitude, 
then the limiting solution as e —> is given by the unperturbed solution of the parabolic 
equation where q has been set to 0. The analysis of the random fluctuations beyond the 
unperturbed solution were addressed in e.g. [H [7]. 

2 Stratonovich integrals and Duhamel solutions 

The analysis of (jH]) requires that we deflne the multi-parameter Stratonovich integral 
used in the construction of a solution to the SPDE. The construction of Stratonovich 
integrals and their relationships to Ito integrals is well-studied. The refer the reader to 
e.g. [6l [TOl [121 [ini [19] . The construction that we use below closely follows the functional 
setting presented in [13] . The convergence of processes to multiple Stratonovich integrals 
may be found in e.g. [31 [1]. 

Let f{xi, . . . , Xn) be a function of n variables in M'^. We want to deflne the iterated 
Stratonovich integral Tn{,f)- Let us flrst assume that / separates as a product of n 
functions deflned on M*^, i.e., /(xi, . . . , Xn) = Y\jt=i fk{xk)- Then we deflne 

n n 

1n{\{hM) = X{Ti{h{xk)). (7) 

k=l k=l 

where Xi(/) = Jj^^ f{x)dW{x) is the usual multi-parameter Ito integral. It then remains 
to extend this deflnition of the integral to more general functions f{x). 
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We define the symmetrized function 



{f,g)r, = E\ I f{x)dW{xi)...dW{xn) g{x)dW{x,)...dW{xr. 

f{x)g{x')E{dW{x^) . ..dWMdW{Xn+l) • ..dW{x2n)}, 



sG6„ 

wfiere tfie sum is taken over tlie n\ permutations of the variables We then 

define Tn{f) = ^n{fs) and thus now consider functions that are symmetric in their 
arguments. 

For the rest of the paper, we write Stratonovich integrals using the notation dW 
rather than o dW. For the Ito convention of integration, we use the notation 6W. Let 
/ and g be two functions of n variables. We formally define the inner product 



(9) 



since the latter has to hold for functions defined as in ([7]). Here, x' = (x„+i, . . . ,X2n)- 
We need to expand the moment of order 2n of Gaussian random variables. The moment 
is defined as follows: 

2n 

E{ Y\ dW{xk)} = XI n ~ xi(^k))dxkdxi(k)- (10) 

k=i pe?!fceAo(p) 

Here, p runs over all possible pairings of 2n variables. There are 

<=-d('P)=c. = JH^i^ = i|ll = (2n-l)!! (11) 

such pairings. Each pairing is defined by a map / = l{p) constructed as follows. The 
domain of definition of / is the subset Aq = Aq{p) of {1, . . . , 2n} and the image of / is 
Bq = Bq{p) = 1{Aq) defined as the complement of in {1, . . . , 2n}. The cardinality of 
^0 and -Bo is thus n and there are c„ choices of the function / such that l{k) > k+1. The 
formula ffTOj) thus generalizes the case n = 1, where 'K{dW{x)dW{y)} = S{x — y)dxdy. 

We extend by density the iterated Stratonovich integral defined in ([7]) to the Banach 
space Bn of functions / that are bounded for the norm 

( / \f f\{Xl, ■ ■ ■ ,X2n) TT 6{Xk - Xil^k))dXkdXi(^k))^ ■ (12) 

The above Banach space may be constructed as the completion of smooth functions 
with compact support for the above norm [18]. Since the sum of product of functions 
of one (i— dimensional variable are dense in the space of continuous functions, they are 
dense in the above Banach space and the Stratonovich integral is thus defined for such 
integrands f{x). A more explicit expression may be obtained for the above norm for 
functions f{x) that are symmetric in their arguments. Since we do not use the explicit 
expression in this paper, we shall not derive it explicitly. We note however that 

= E{X„H_„(|/®/|)}, (13) 
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since J„(/)J„(/) = J2n(/ ® /)• 

Note that the above space is a Banach subspace of the Hilbert space of square 
integrable functions since the norm of / appears for the pairing Hfe ^(^fc ~ Xk+n)- 
Note also that the above space is dense in L'^(W^'^) for its natural norm. Indeed, let 
/ be a square integrable function. We can construct a sequence of functions /'^ that 
vanish on a set of measure in the vicinity of the sets of measure where the 
distributions S{xk — xi), 1 < k,l < n, are supported and equal to / outside of this 
set. For such functions, we verify that ||/'^||n is the L^(]R"'^) norm of Moreover, 

converges to / as /c ^ oo as an application of the dominated Lebesgue convergence 
theorem so that i3„ is dense in L^(]R"'^). Note finally that the above expression still 
defines a norm for functions that are not necessarily symmetric in their arguments. 
This norm applied to non-symmetric functions is not optimal as far as the definition 
of iterated Stratonovich integrals are concerned since many cancellations may happen 
by symmetrization ([H]). However, the above norm is sufficient in the construction of a 
Duhamel expansion solution to the SPDE. 

Duhamel solution. Let us define formally the integral 

Hu{t,x)=a[ [ G{t- s,x]y)u{s,y)dW{y)ds, (14) 

Jo JR'i 

where we recall that dW means an integral in the Stratonovich sense. The Duhamel 
solution is defined formally as 

oo 

= Un+l{t,x) =nUn{t,x), Mq (t, x) = e"*-^^'''^) [Mq (x)] , (15) 

n=0 

where Uq is the initial conditions of the stochastic equation, which we assume is inte- 
grable. The above solution is thus defined formally as a sum of iterated Stratonovich 
integrals Un(t, x) = In{ fn{t, x,-)). 

The main result of this section is the following. 

Theorem 1 Let u{t, x) be the function defined in f[T^ . The iterated integrals Un{t, x) = 
^nifn{t,x, ■)) are defined in L'^(W^; Bn) uniformly in time t G (0,T) for all T > 
and n > 1. When the initial condition uo{x) G L-'^(]R'^) fl L^(M'^), then u{t,x) is a 
mild solution to the SPDE in L^{M.'^ x Q) uniformly in time t G (0,T) for all T > 0. 
When Uq{x) G L^(R'^), then the deterministic component UQ{t,x) in u{t,x) satisfies 
t2^Uo{t,x) G L^(M'^) uniformly in time. 

Proof. The norm of u{t, x) is defined by 

/ E{u'^{t,x)}dx = V / E{In+mifnit,X,-)0 fm{t,X,-))}dx 

n,m>0 

< J2 f ^{In+m{\fn{t,X,-)® fm{t,X,-)\)}dx. 

n,m>0 

We now prove that the latter is bounded uniformly in time on compact intervals. The 
proof shows that fn{t,-) is also uniformly bounded in L^(M'^;i3„) so that the iterated 
integrals Un{t,x) are indeed well defined. 
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n-l .tk 



n 



m—1 



Note that n + m = 2n for otherwise the above integral vanishes. Then, using the 
notation = sq = t, we have 

In,m{t) = / E{In+m{\fn{t,x) fm{t,x)\)}dx = 

n—1 „ n 

W\G\{tk-tk+i,Xk;Xk+i) I G{tn-,Xn]i)uQ{C}di\Wdtk 

to ' tl 

n/ / ni^l'^^^~*'+i'^''^'+i) / G{sm,yni]QuQ{Cidc\Wdsi 

n m 

d{xo - x)6{yo - x)a"+™E{ JJ dW{xk) JJ diy(yz)}rfa;. 

fc=i 1=1 

Using the fact that 2ab < + b"^ with a and b the Green's functions involving x and 
the fact that r m J (^^(t, x, is uniformly bounded, we bound the integral in x by a 
constant. Let us define (p{s) = |t — s|~m . As a consequence, we obtain that 



n-l 





m—1 

n 



k=l 




JR<*("-1) 
m—1 



k=l 



Y\_\G\{tk -tk+i,Xk;xk+i) / G{tn,Xn;Ouo{OdnY\_dt 



k=l 



JE''(™^ 
n m 



III, — J_ n III' 

Y\ \G\{si- si+i,yi]yi+i) / (^(s^, ?/m; C)%(C)t^C TT 

m 

a"+-E{ J] dW{xk) n 
fc=i «=i 

Here a < b means that a < Cb for some constant C > 0. 

Let us re-label Xn+i = yi and tn+i = s/ for 1 < / < m. We also define x 
{xq, . . . , Xn+m+i)- Then we find that 

2n 



< / Hr,,m{t, x) e| ft dW{Xk) 

/■t n—1 m—1 „ 

H^,m{t,^)= / 0(ti)n / / n / 

-'0 -'0 -'0 Jo 



n—1 „f, „f m—1 „f , , , , n+m— 1 

, Xk , Xk-\-l ) 

k=l,kj^n 

n+m 



G{tn,Xn;Ouo{C)di I G{tn+m-:Xn+m-X)Uo{C)dQW dtk- 

' Li 

We now recall the pairings introduced in (|TOl) and replace n by n there. Let us introduce 
the notation 



Xk+i k n,n + m 
yk= ■{ ^ k = n 

C k = n + m 

so that Hn,m{t,^ is bounded by 



n-l „tk ftm-l .t„+i 2n 

<^(^i)n / / n / ni^K^ 



k n,n + m 
k = n,n + m, 



]R2d Jo 



fc=l 



Jo 



«=1 



n+m 



Cfc - Tk,Xk;yk)\uo 



(OII«o(C)l«n^4- 



fc=i 



fc=i 



Now, we have for each pairing p G 

2n 

W\G\{tk - Tk,Xk,yk) = Yl \G\itk - Tk,Xk;yk)\G\{ti(^k) - n{k),xi(^k);yi{k)), 

k=l keAo 

and as a consequence, using the delta functions appearing in (fTIl . 

n— 1 ^f, m—1 



in,ut) < / / '^(^i) n / / n / i^o(oi i^o(c)i 

„ n+m 

"Tk ) Xk ] yk)\G\{ti(^k) - Ti(^k),Xk;yi(k))dxkj d^dC\\dtk (16) 



fce^o k=i 

s E / ^(*o n / / n / n n^oinn*- 

for some positive constant C in which we absorb . On the second hne above, the yi(k) 
are evaluated at xi{k) = Xk- The function k ^ t{k) for k & Aq is at the moment an 
arbitrary function such that t{k) = k or i{k) = l{k). The last line is obtained iteratively 
in increasing values of k in Aq by using that one of the Green's function is integrable 
in Xk uniformly in the other variables and that the other Green's function is bounded 
independent of the spatial variables by a constant times the time variable to the power 
—a with a := —. We have used here assumption (jl]). It then remains to integrate in 
the variables ^ and C, and we use the initial condition uq{x) for this. 

Let us now choose the map t{k). It is constructed as follows. When both k and l{k) 
belong to {1, . . . , n} or both belong to {ri + 1, n + m}, then we set t{k) = k. When 
/c G {1, . . . , n} and l{k) G {n + 1, n + m} (i.e., when there is a crossing from the n first 
variables to the m last variables), then we choose i{k) = k for half of these crossings and 
t{k) = l{k) for the other half. When the number of crossings is odd, the last crossing is 
chosen with t{k) = k. 

Let us define Al = i{Ao) n {0, . . . ,n} and Al = t{Ao)\Al. Let uq = no{p) be the 
number of elements in A^ and mo = mo(p) be the number of elements in Aq such that 
Uq + rriQ = n. Let p = p(p) be the number of crossings in p. Then, by construction of 
m, we have 



n — p 
^0 = —i; \- 



p + l 



m — p 
"^0 = — t; — 



(17) 



where = 2±i if p is odd and | if p is even, with [^^] + [|] = p. Thus, riQ is bounded 
by ^ and rriQ by y. 



We thus obtain that 

n-l 



m 



peqj A;=o 

m—l 



[n 



1=0 n+l&Al ^'+^' 1=1 



with the convention that = Sq = t, = and Sm+i = 0. It remains to estimate the 
time integrals, which are very small, and sum over a very large number of them. It turns 
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out that these integrals admit exphcit expressions. The construction of the mapping 
t{k) ensures that the number of singular terms of the form r~" is not too large in the 
integrals over the t and the s variables. 

Let CKfc for < k < n he defined such that = a, ak = a for k ^ and ctfc = 
otherwise. Still with the convention that tn+i = 0, we thus want to estimate 

fc=0 -^^ k=0 ^ '^^'-^ k=l 
The integrals are calculated as follows. Let us consider the last integral: 

-du, 



1 - 

where we define /3„ = a„, and = Pm+i + ctm for < m < n — 1. The latter integral 
is thus given by 

where B{x,y) = ^^^^^y-^ is the Beta function and T(x) the Gamma function equal to 
(x — 1)! for a: G N*. The integration in t„„2 then yields 

/ 77 ~4 ^^~^'^^"-2 = ^n-2'^B{2 - (3n-i, 1 - an-2)- 

Jo [tn-2 - tn-l)""-2 

By induction, we thus obtain that 

Since + = /?fc, we obtain by telescopic cancellations that 

n-l 



^ ^^-/^o r(i - Pn) TT r(i _ 



k=0 

Then with our explicit choices for the coefficients ak above, we find that Po = (no + l)a 
so that 

J ^ ^n-(no+l)a r(l - an)r"''(l - a) 

"~ ° r(n + l-(no + !)«)■ 

For a fixed p, we see that the contribution of the time integrals in In,m{i) is bounded 
by a constant (since r(l — a) is bounded as a < 1) times 

r"(l-a) ^ r"(l-a) 



r(n + 1 - (no + l)Q!)r(m + 1 - rriQa) " r((n + 1)(1 - f ) - a)r((m + 1)(1 - f )) 

based on the values of riQ and mo. Using Stirling's formula T{z) ~ (^)^{^Y so that 
r(2;) is bounded from below by {^Y for C < e, we find that the latter term is bounded 
by 

^n(l-f)^m(l-f)' (20) 



for some positive constant C. The latter bound holds for each p G Using the Stirling 
formula again, we observe that the number of graphs in ^ is bounded by (— As a 
consequence, we have 

Using the concavity of the log function, we have 



vJ'm^ > ( ■ ) > 

n + m 



so that 



As a consequence, we have 



T ^ J . ,(n+m){l 2.) a„^„^ 1 (')0\ 

J-n,m ^ >Jn,m\>-) ■— 1^0 znci-ryV ^ ' 

The bound with n = m shows that for n > 1, M„(t, x) belongs to L^(M'^; Bn) uniformly 
in time on compact intervals since 2(1 — a) > 0. Now the deterministic component 
Uo{t, x) is in L^(]R'^) uniformly in time when Uq{x) G L'^{W^) while tiuoit, x) is in //^(M*^) 
uniformly in time when Uq{x) G L^(R'^). Upon summing the above bound over n and 
m, we indeed deduce that u{t,x) belongs to L^{Q x Mf^) uniformly in time on compact 
intervals when uq G L^(]R'^). 

The above uniform convergence shows that T-tu{t,x) is well defined in L^(f2 x M.'^) 
uniformly in time. Moreover, we verify that Huit, x) = ^„>i Un(t, x) = u(t, x)—Uo{t, x). 
This shows that u{t, x) is a mild solution of the stochastic partial differential equation 
and concludes the proof of the theorem. □ 

3 Uniqueness of the SPDE solution 

Let us assume that two solutions exist in a linear vector space DJl. Then their difference, 
which we call u, solves the equation 

u = Tiu = H^u, 

for all p > 0. The space DJl is defined so that Ti^u is well-defined and is constructed as 
follows. 

We construct w G 971 as a sum of iterated Stratonovich integrals 



u{t,x) = ^ J„(/„(t,x, ■))• 



m>0 



Because the iterated Stratonovich integrals do not form an orthogonal basis of random 
variables in L^(f2), the above sum is formal and needs to be defined carefully. We need 
to ensure that the sum converges in an appropriate sense and that 9Jt is closed under 
the application of H. 
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One way to do so is to construct u{t, x) using the classical Wiener-Ito chaos expansion 



U{t,x) = ^ Imigm{t,X, ■)), 

where is the iterated Ito integral, and to show that the above series is well defined. 
We then also impose that the chaos expansion of Ti^u is also well-defined. 

We first need a calculus to change variables from a definition in terms of iterated 
Stratonovich integrals to one in terms of iterated Ito integrals. This is done by using 
the Hu-Meyer formulas. We re-derive this expression as follows. We denote by 6W an 
Ito integral and by dW a Stratonovich integral. We project Stratonovich integrals onto 
the orthogonal basis of Ito integrals as follows 

E{In{ fn)Im{<Pm)} = lE{/„ (0^) } = ml gm4>mdx, 

where 0^, is a test function. We find that K{Xn{fn)Im{4'm)} is equal to 

fn{Xi, . . . , Xn)(j)m{yu ym)E{dW{Xi) . . . dW{Xn)SW{yi) . . . SW{ym)}. 

[7i~\-m)d 

The moment of product of Gaussian variables is handled as in (fTOl) with the exception 
that K{5W {yk)SW (yi)} = for k ^ I hj renormalization of the Ito-Skorohod integral. 
The functions /„ and (pm are symmetric in their arguments (i.e., invariant by permuta- 
tion of its variables). We observe that the variables y need be paired with m variables 
X. There are ways of pairing the y variables. There remain n — m = 2k variables 
that need be paired, for a possible number of pairings equal to 

(2A;- 1)! 



(A; - l)!2^-i' 

The above term is thus given by 

m + 2k\ {2k -1] 
m 



{k - 1)!2^"^ J ij f'^+^l<^(y^^---yrni^li^l^---^k,Xk)YldXl^(prn{yi,---ym)Yldy 

This shows that Qm is given by 

, , (m + 2k)l f , , ^2 A , 

gm[xi, ■ ■ . = J fm+2k{xi, . . . , • • • , l/fc ) [[dyk- 



p=i 



Here ?/®^ = {y,y)- The coefficients Qm are therefore obtained by integrating n — m 
factors pairwise in the coefficients /„. This allows us to write the iterated Stratonovich 
integral as a sum of iterated Ito integrals as follows: 

f-l 
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This is the Hu-Meyer formula. More interesting for us is the reverse change of coordi- 
nates. Let us define formally 

/ = ^^nifn) = ^ Im{9m)- 

n>0 m>0 

Then we find that 

^ (m + 2k)\ f , m\ 1 

9m{x) = > — ^ / fm+2k{x,y"' )dy. 

The square integrability of the coefficients gm is a necessary condition for the random 
variables / to be square integrable, and more generally, to be in larger spaces of distri- 
butions [8]. The above formula provide the type of constraints we need to impose on 
the traces of the coefficients /„. For square integrable variables, we consider the normed 
vector space 9Jt/ of random variables 

n>0 

where the coefficients {/n} are bounded for the norm 

Note that the above defines a norm as the triangle inequality is clearly satisfied and for 
= 0, we find that the norm of each has to vanish, so that = for all m 
when the norm vanishes. Note also that 971/ is a dense subset of LP'{VL) as any square 
integrable function may be approximated by a function which vanishes in a set of 
Lebesgue measure at most k^^ in the vicinity of the measure set of diagonals given by 
the support of the distributions 6{xk — xi). For such functions, we verify that = g^ so 
that the Ito and Stratonovich iterated integrals agree. We also have that g^ converges 
to gm by density. Since every square integrable random variable may be approximated 
by a finite number of terms in the chaos expansion, this concludes our proof that SOT/ is 
dense in L?{VL) equipped with its natural metric. 

Let us now move to the analysis of the stochastic integral Ti. It turns out that 971/ 
is not stable under Ti nor is it in any natural generalization of 971/. Let us define 

u{t,x) = ^ J„(/„(t,a:, ■)), Huit.x) = ^In{{nf)nit, x, ■)) 

n>0 n>0 

We then observe that 

nfn+i{t,x,y) = as[ / G{t - s,x,yi)fn{s,y)ds], 
Jo 

where 5 is the symmetrization with respect to the d{n + 1)— dimensional y variables. 
Let us consider 7i^/„+2, which depends only on /„. Let Cm,k = ^mtm^' coefficient 
that appears in the definition of gm- Then, for Ti? fn+2 relative to /„, the coefficients 
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indexed by k are essentially replaced by coefficients indexed by A; + 1. Since Cm,k+i is not 
bounded by a multiple of Cm,k uniformly, the integral operator Ti^ cannot be bounded in 
9Jl/. The reason why solutions to the stochastic equation may still be found is because 
the integrations in time after n iterations of the integral Ti provide a factor inversely 
proportional to n\. This factors allows us to stabilize the growth in the traces that 
appears by going from Cm,k to Cm,k+i- Uniqueness of the solution may thus only be 
obtained in a space where the factor n! appears, at least implicitly. 

A suitable functional space is constructed as follows. Let Qm be the chaos expansion 
coefficients associated to the coefficients and gm,p the chaos expansion coefficients 
associated to the coefficients |?i^/n|- 

Then we impose that the coefficients {/n} be bounded for the norm 

/I 
g"^ {t,x,y)dxdy] < oo, (23) 

where Cp is an increasing series such that Cp — > oo as p — > oo. Here T is a fixed 
(arbitrary) positive time. We denote hj DJt = OJt(T) the normed vector space of random 
fields u{t, x) for which the decomposition in iterated Stratonovich integrals satisfies the 
above constraint. 

We are now ready to state the main result of this section. 

Theorem 2 Let T > be an arbitrary time anduo{x) G L^(]R'^) nL^(M'^). The solution 
constructed in TheoremU\is the unique mild solution to the stochastic partial differential 
equation ^ in the space 971 = 9Jl(T) . 

Proof. Let us first prove uniqueness in 971. We have u = IHPu for all p > 0. This 
implies that gmit, x, ■) = gm,p{t, x-). The latter converges to in the sense as p ^ oo. 
This implies that gm{t,x, ■) uniformly vanishes for all m so that u = 0. 

Let now iL{t, x) be given by the following Duhamel expansion 

u{t,x) = 'Y^Un{t,x), Un+l{t,x) = T-CUn{t, x) = H^'^'^Uo^t, x) , Uo{t,x) = C"^^ Uq{x) . 

We thus verify that 

H''u{t,x) = y^/Un{t,x). 

n>k 

This shows that the norm of 7i^u(t, ■) is bounded by the sum of the coefficients 
In,m{'t) for n,m > k. This sum clearly converges to as A; ^ oo. Call this sum c^^. We 
recall that 

E{u'^{t,x)}dx = Ym\ j gl^{t,x,y)dxdy = E J ^Xn+mifn® fni){t,x)dx. 



m>0 



Then we find that 



Cfc j gl^^kii, X, y)dxdy < Ck\\u{t, ■)\\l2 
This shows that u belongs to 971. □ 



(f^xRd) — ^■ 
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The same theory holds when the supremum in m is replaced by a sum with weight ml 
so that 971 becomes a subspace of L^. In some sense, the subspace created above is the 
smallest we can consider that is stable under application of 7i. When Uq{x) G L^(]R'^) not 
necessarily in L^(M"'), then the deterministic component UQ{t,x) is not square integrable 
uniformly in time. The space 971 may then be replaced by a different space where Cp in 
(123!) is replaced by t^Cp. 



4 Convergence result 

Let us now come back to the solution of the equation with random coefficients ([1]). The 
theory of existence for such an equation is very similar to that for the stochastic limit. 
We define formally the integral 

neu{t,x)= / G{t- s,x]y)u{s,y)qs{y)dyds, (24) 

Jo JR'i 

where we have defined ^^(y) = e~^q{^). The Duhamel solution is defined formally as 

oo 

Ue{t,x) = '^Un,£{t,x), Un+l,£{t, x) = HeUn,e{t, x) , Uo(t, x) = e"*-^^'^''^) [Mo(a;)] , (25) 



n=0 



where uq is the initial conditions of the stochastic equation, which we assume is inte- 
grable. We have the first result: 

Theorem 3 The function u^{t,x) defined in (1251) solves 

u,{t,x) = HMt.x) + e-*^(^'^)K(x)], (26) 
and is in L^(]R'^ x Q) uniformly in time t G (0, T) for all T > 0. 

Proof. The proof goes along the same lines as that of Theorem [H The norm of 
Ue(t,x) is defined by 

/ E{ul{t,x)}dx = ^ E{Un,£{t,x)Um,e{t,x)}dx 

< / ^{\Un,e\{t,x)\Um,e\{t,x)}dx < I^^n,e{^), 



n,m>0 ' 



where 



In,m,s{t)=l Wl I W\G\{tk-tk+l,Xk;Xk+l) 1 G{tn,Xn;0'^o{C)di\Wdtk 



Wl I TTl'^l(^'^^'+i'^''^'+i) / G{s^,ym]Ouo{Qdc\Wds 
to J^""- 7=0 ' -^^^ ' 7=1 

n m 

5{xo - x)S{yo - x)E{Y\_ qe{xk)dxk W qe{yi)dyi}dx. 



fe=i 1=1 
13 



Following the proof of Theorem [T], we obtain 



The statement (fTOj) now becomes 

2n 

E[l[q,{xk)dxk] = Y, n (27) 

k=i peq3fceAo(p) 

where we recall that R{x) = K{q{0)q{x)} is the correlation function of the Gaussian 
field q. This yields 

(• ft ftk ft ftn + l f 

in,mAt) < E / / '^(^i) n / / n / i^o(^)i / n 

~:zJK2dJ(, f- ,Jo Jo Jo Jr"'' ^^-^ 



Xk — Xi{k) , 



G\{tk - Tk,Xk]yk)\G\{ti{k) - ri(k),xi(^k)]yi(k))e R{ )\dxkdxi(^k)) d^dCY[dtk 



n+m 



k=l 



For each k & Aq considered iteratively with increasing order, the term between paren- 
theses is bounded by the norm of the Green's function integrated in times the 
integral of the correlation function in the variable with {t{k),t'{k)) = {k,l{k)), 

which gives a contribution thanks to the definition ([3]), times the L°° norm of the 
Green's function in the variable Xf'(fc). Using (jlj) and the integrability of the correlation 
function R{x), this shows that 

ft '^-i ftk ft"^-'^ ft„+i "+"1 

w(t) < E / '^(^i) n / / n / n ^ wuowi.wdtk, 

peqj-^o fc=i>^o -^0 «=Wo k&Ao \Hk) - ^{k)) ^ k=i 

as in the proof of Theorem [H The rest of the proof is therefore as in Theorem [1] and 
shows that each is well defined in L'^{W^ x ^2) uniformly in time and that the 

series defining u{t, x) converges uniformly in the same sense. □ 



Mollification and convergence result. We now have defined a sequence of solu- 
tions Us{t,x) and a limiting solution u(t,x). When q^ and the white noise W used in 
the construction of u(t, x) are independent, then the best we can hope for is that 
converges in distribution to u. The convergence is in fact much stronger (path-wise) by 
constructing qedx as a moUifier of dW. Let -R(0 be the power spectrum of q, which is 
defined as the Fourier transform of R{x). By Bochner's theorem, the power spectrum 
is non- negative and we may define p{^) = {R{^))2. Let p{x) be the inverse Fourier 
transform of p. We may then define 

q{x) = I p{x-y)dW{y), (28) 

and obtain a stationary Gaussian process q{x). This process is mean-zero and its cor- 
relation function is given by 

R{x) = p{x - y)p{y)dy = R{x), 

jRd 
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by inverse Fourier transform of a product. As a consequence, q{x) and q{x) have the same 
law since they are mean zero and their correlation functions agree. The corresponding 
Duhamel solutions and £te also have the same law by inspection. It thus obviously 
remains to understand the limiting law of to obtain that of u^. It turns out that 
may be interpreted as a moUifier of u{t, x), the solution constructed in Theorem [U and 
as such converges strongly to its limit. 

In addition to the assumptions on the Green's function in (jlj) and ([5]), we also assume 
that p{x) G L^(M'^). Then we have 

Theorem 4 Let Us{t,x) be the solution constructed in Theorem\^ and u{t,x) the solu- 
tion constructed in TheoremUi Then we have that Us{t,x) converges in distribution to 
u{t,x) as £ — ^ 0. More precisely, let Ue{t,x) be the Duhamel solution corresponding to 
the random potential q in f l28l) . Then we have that 

-M(t)||i2(iKdxQ) 0, £-^0, (29) 

uniformly in time over compact intervals. 

Proof. Let us drop the upper" to simplify notation. We have 

61,{t) = [ E{{u{t)-u,{t)f}dx = y^6I,,n,M 



Following the proofs of Theorems [T] and [3l we observe that 

n— 1 „(j, „ n— 1 „ n 



TT / / TT'^(^fc / G{tn,Xn;Ouo{Od^Y\dtk 

k=J^ ■J^'-k=o k=l 

m—1 „ m—1 „ m 

TT / / TT '^(■^^ ~" / G{s^,ym;C)uo{OdCY\dsi6{xo-x) 

to Jo J^'"" 7=0 -J^' 7=1 

n n mm 

6{yo - 2;)e| JJ adW{xk) - JJ qe{xk)dxk^ JJ (TdW{yk) - JJ qe{yi)dy?j }c?x. 



k=l k=l 1=1 



Here, we have again that to = sq = t. The integration in x is handled as in the proof of 
Theorem [T] so that 



\5Ie,n,m{t)\ < 

-t n-1 ^ n-1 



^(^i)TT/ / T] G{tk- tk+i,Xk;xk+i) G{tr„Xn]Ouo{OdCY{ dtk 

7=1 Jo J^''-''k=i 7=1 

„t m—1 „ m—1 „ m 

/ TT / / W G{si- si^i,yuyi+i) I C)%(C)c^C TT ^'^^ 

^0 iJ{ Jo jRd(m-i) J^d 

n n m. m 

e{ ( H adW{xk) - n Qs{xk)dxk) ( n ^^^W^(yfc) - TT ^e{yi)dyi 



k=l k=l 1=1 1=1 
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The main difference with respect to previous proofs is that we cannot bound the Green's 
functions by their absolute values just yet. The moment of Gaussian variables is handled 
as follows. We recast it as 



Y\ adW{xk) - Y\ qe{xk)dxk] Y\ (le{yi)dy, 



fc=l k=l 1=1 



plus a second contribution that is handled similarly. We denote by 51] ^^{t) the corre- 
sponding contribution in 5Ie,n,m{f) and by = SIe,n,m{t) — 5ll.^^^{t). The above 
contribution is recast as 

n g— 1 n+m 

^^adW{xk)(adW{Xq) - qe{xq)dx^ JJ qe{xp)dxp, (30) 

5=1 p=l p=q+l 

where we have defined Xn+i = yi ioT 1 < I < m. We have therefore n (or more precisely 
n A m by decomposing the product over m variables when m < n) terms of the form 

g— 1 n+m 2n 

Ie| Y\_^dW{xk)(^adW{xq) - qe{xq)dx^ J]^ qe{xp)dxp^ := e| J]^ afc,£((ixfc)|, 

p=l p=q+l k=l 

where each measure ak^e{dxk) is Gaussian. Then, f fTOj) is replaced in this context by 

2n 

Ie{ JJafc,e(c^a;fc)| = X] 11 '^{o'k,e{d^k)ai{k),e{d^m)] 

k=i pGq3fceAo(p) (31) 

:= ^ JJ K,kixk - Xl(^k))dxkdxl(^k)■ 
pe<^keAo{p) 

The functions h^^ki^k — xi(^k)) come in five different forms according as 

E{dW{x)dW{y)} = a^{x-y)dxdy 

1 /X — y 



E{dW{x)q,{y)dy} = a-^p(^^—^yxdy 
E{q,{x)dxq,{y)dy} = —^R[-^—Jdxdy (32) 



E{{dW{x)-qe{x)dx)dW{y)} = {a^5{x - y) - (y^p{^-j^y^dxdy 
E{{dW{x) - q,{x)dx)qMdy] = (^^p(^) " J,R{^))dxdy. 



At this point, we have obtained that 



4>{h)W I / WG{tk-tk+i,Xk]Xk+i) I G{tn,Xn]CjUo{OdiWdtk 

TT / / TT ^(^' ~ / G{sm,ym;Ouo{OdcY\^^i 

T=i ti ti 

Y\_ K,k{xk - xi(k))dxkdxi(^k) 

fc6Ao(p) 
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Using the notation as in the proof of Theorem [H we obtain that 

winAt)\<Y. \ / / '^(^on / / n / «o(eH(c) / n 

n+m 



G{tk - n,Xk]yk)G{ti(^k) - n(k),xi(ky,yi{k))hs,k{xk - xii^k))dxkdxi{kyjd^dC, J]^ dt 



A: I 

fc=l 



<E / / '^(^on / / n / ivoiMoi / n 

n+m 

G{tk - Tk,Xk]yk)G{ti(^k) - n(k),xi(^ky,yi{k))he,k{xk - xi(k))dxkdxi^k))\d^dC JJ dtk. 



k=l 



It remains to handle the multiple integral between absolute values. For k G ^o(p) 
for which h^ ^ is of the form given in the last two lines of fl32l) . we observe that the 
corresponding term between parentheses in the above expression is of the form 



G{s, x; OG{t, y; OMx - y)dxdy 

= j^^ G{s, X, C) ( ^9{^) (G(r, x; - G(r, y; 0)dy)dx (33) 

= / G{s,xX)( g{y){G{T,x;^) -G{t,x + ey;0)dy)dx, 

where the function g{x) is given by 

g{x) = ±ap{x), or g{x) = ±{R{x) — ap{x)). (34) 

This is because p averages to a while R averages to a^. Let kQ be the index for which 
h^^ko is in the form of a difference as above. This yields, with g = g[kQ] as above, 

m,n,mit)\<Yi / / <^(^i)n / H I K(^)ii^o(c)i / n 



G{tk — Tk, Xk] yk)G{ti(^k) — Ti{k),xi(k); yi(k))he,k{xk — xi(k)] 



dxi-dx 



n+m 



\G{tii^ko) -Ti(^ko)-,Xt{koyyi(ko))\\9{xt'(ko))\d^dC JJ dtk 



{G{-, -, ■) - G(-, ■ + exiif^ko), ■))(^r(fco) - Tr(feo)' ^e(fco); Vviko)) 



kaxi(k) 



dXkodXK^ko)- 



The above term is now handled as in the proof of Theorem [H For k ^ ko, the bounds 
are obtained as before because p and R are integrable functions by hypothesis. The 
Green's function \G{ti(^ko) ~ ''"e(fco)' ^t(A:o)j l/e(A:o))l is bounded by a constant times |tt(A:o) — 
'^t{ko)\~°'- The integration dxk^^dxK^ko) = dxi(^ko)dxii(^ko) then yields a contribution bounded 
by |tr(feo) - Trcfco)!"'^ times 

Me= sup f \g{y)\\G{T,x-i)-G{T,x + ey-0\dxdy. (35) 
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The presence of the factor 7 is necessary in order for Mg to converge to as e ^ 0. As 
a consequence, as in the derivation of f|T6l) . we obtain that 

< E / n / / n / ^ — ^^i^^ n 7. ^n^^'^- 



The factor 2n comes from twice the summed contributions in fl30|) . The presence of the 
factor 7 increases the time integrals as follows. Assume that < n for concreteness; 
the case A;o > n + 1 is handled similarly. Then /3o in the proof of Theorem [T] should be 
replaced by (3q + 7. This does not significantly modify the analysis of the F functions 
and the contribution of each graph is still bounded by a term of the form (!20|) . The 
behavior in time, however, is modified by the presence of the contribution 7 and we find 
that 



„t(l-a)^f(l-a)- 

The above bound is of interest for n + m > 2 since the case n = m = corresponds 
to the ballistic component UQ{t,x), which is the same for Ue(t,x) and u{t,x) so that 
^h,o,o = 0. By choosing 7 = 2(1 — a) > 0, we observe that 2(1 — |) — a — 7 > for 
n + m > 2 so that |(^/e,n,m(OI is bounded uniformly in time. The new factor n may be 
absorbed into so that after summation over n and m, we get 

\K{t)-u{t)\\l^^,^^^<CMe. (36) 

By assumption (Q, the integrand in (l35ll converges point-wise to and an application 
of the dominated Lebesgue convergence theorem shows that — > 0. This concludes 
the proof of the convergence result. □ 



A continuity lemma. We conclude this paper by showing that the constraints (jlj) 
and ([5]) imposed on the Green's functions of the unperturbed problem throughout the 
paper are satisfied for a natural class of parabolic operators. 

Lemma 4.1 Let G{t,x) be defined as the Fourier transform o/e"*'^'", i.e., the Green's 
function of the operator dt + (—A)" for m> d. Then the conditions in (jlj) and ^ are 
satisfied. Moreover, when m is an even number, then in (15^ satisfies the bound 

M, /3 = 2(m-t/) Al. 

Proof. By scaling, we find that G(t,x) = t~'^G{l,t~^x). Since l^pe"'^'"" isintegrable 
for all p, we obtain that G{l,x) belongs to C°°(M'^). Since G{l,x) is bounded, then so 
is t'^G{t.,x) uniformly in t and x. 

By the above scaling, G{t,x) belongs to L^(]R'^) uniformly in time if and only if 
^(l,^) does. When m in an even integer, then e~'^'™ belongs to 5(M'^), the space of 
Schwartz functions, so that G{l,x) G iS(M'^) as well. It is therefore integrable and has 
an integrable gradient. When m is not an even integer, we have 

e-i«r-i = 5^^iei'=-. 

k=l 
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The Fourier transform of the homogeneous function |.^|'^'" is given by [20 
c(/c)|xr^"^-^ c{k) = Cd2 



-km-d ^ofcmr(i(^m + rf)) 



where Cd is a normahzation constant independent of k. The Fourier transform of e~'^'" 
may then be written as a constant times plus a smoother contribution that 

converges faster to (for instance because it belongs to some H'{R'^) with s > ^ + k 
sufficiently large so that k derivatives of this contribution are integrable). It is therefore 
integrable for m > 0. The bound follows from the and L°° bounds. 
We obtain by scaling and from definition of G{t, x) that 

P \G{t,x) -G{t,x + ey)\dx = f' \G{l,x) - G{l,x + r^ey)\dx. 

The above derivation shows that the gradient of G is also integrable for m > 1 so we 
may bound the above quantity by P{1 At~'^e\y\). Now, 

snp{r Ar~^e\y\) < ie\y\r^y e\y\ = ie\y\f^'^~''^Ve\y\, 

t<T 

according as 7m < 1 or 7m > 1. We thus obtain ([5]) by sending ey 0. When g{y) is 
sufficiently regular, then we obtain the more precise bound 

M^<e'i--^) ! \g{y)\\y\'^-~''^dyye[ \g{y)\\y\dy, 

provided that the latter integrals are well-defined. □ 
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